Spatial distribution of atoms in gas-covered Pd-X 
nanoparticles (X= Ag, Cu, Ni, Pt) 

Mahesh Menon and Badal C. Khanra* 
Condensed Matter Physics Group, Saha Institute of Nuclear Physics 
1/AF, Bidhannagar, Calcutta- 700 064, India 

Abstract 

A Monte-carlo (MC) simulation procedure has been developed where the pair 
bond energies are allowed to take into account the various coordination numbers of 
surface atoms and the presence of adsorbates. The pair bond energies are calculated 
from partial bond energies of atoms which, in turn, are calculated from modified 
tight binding model in the second moment approximation. The model has been 
applied to study the role of adsorption of hydrogen, oxygen, carbon monoxide and 
nitric oxide on the surface composition and surface bond geometry of bimetallic 
Pd-X (X = Ag, Cu, Ni, Pt) nanoparticles having fee cubo-octahedral geometry 
with 201, 586, 1289 and 2406 atoms. The results are compared with the known 
experimental results. Importance of the results in studying reactions on supported 
bimetallic catalysts has been highlighted. 
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I. INTRODUCTION 

Supported bimetallic catalysts are of enormous importance in petroleum and 
other chemical industries. They are used for better selectivity of products in 
parallel reactions. Mostly, these bimetallic catalysts are supported nanostructures, 
which may have specific active sites like the corner sites, edge sites or sites on reg- 
ular crystalline faces. In order to understand the catalytic activity and selectivity of 
a bimetallic nanostructure it is important to know the average surface composition 
and particularly the occupancy of various sites by A or B metal atoms of the A-B 
bimetallic system. Since in catalytic reactions the bimetallic catalyst surfaces are 
partially or totally covered by reacting gases, it is also important to study the sur- 
face composition in presence of reacting gases. It is the aim of the present work to 
study the spatial distribution of atoms for a few Pd-based bimetallic nanoparticles 
like Pd-Ag, Pd-Cu, Pd-Ni and Pd-Pt nanoparticles having dispersion in the range 
of 0.6 - 0.3 (the total number of atoms per particle being in the range of 200 - 
2500). The motivation to study Pd-based systems follows from the fact that Pd is 
a very good catalyst for hydrogenation and oxidation reactions. Amongst the alloy 
systems, Pd-Pt catalysts are of great technological interest because of their use in 
(de)hydrogenation reactions, oxidation reaction, (dehydro)cyclization reaction, hy- 
drogen peroxide production etc,"^"^ besides their use in electrochemical studies. 
These catalysts have the additional advantage of being more resistant to poisoning 
by sulfur and nitrogen compounds than the respective monometallic catalysts. 
Alloying Pd with copper and silver (Pd-Cu, and Pd-Ag systems) is found to greatly 
enhance the selectivity of Pd in the partial hydrogenation of dienes^^ and in the 
hydro-dehalogenation of olefinic compounds.^^ In these systems the Pd surface has 
to be partly covered by the second metal (Cu or Ag) in a carefully controlled way. 
This is achieved in bimetallic nanostructures. Besides, Pd-Cu systems show some 
promise as catalysts for simultaneous oxidation of CO and reduction of NO,^''"^^ 
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the very vital reactions for automotive pollution control. Pd-Ni systems also show 
interesting features in their use as catalysts for 1,3-butadiene hydrogenation reac- 
tion. The PdsNigs system, for example, is found to have greater activity than the 
pure metals. In view of the above interesting features we study in this work the 
atomic distribution in the above four Pd-based nanostructures in presence of differ- 
ent coverage (0 < 6* < 1) of some simple adsorbates like hydrogen, oxygen, carbon 
monoxide (CO) and nitric oxide (NO). The plan of the paper is as follows: In section 
II we describe the formulation of the model. In section III we apply the formalism 
to the four Pd-based bimetallic nanoparticle system mentioned above. Conclusions 
are drawn in section IV. 

II. FORMULATION OF THE MODEL 

The model used in this work is based on the Monte-Carlo (MC) simulation 
procedure adopted by Strohl and King^^ to study the segregation behavior of Pt-lb 
bimetallic nanoparticles. However, we have introduced two modifications to their 
procedure as discussed below: 

First, the pair bond energy required for MC simulation are evaluated taking into 
account the fact that the pair bond energies are dependent on the coordination 
number of the atoms. Thus, the pair bond energy between two nearest neighbor 
atoms, Eij, is related to the site energy (or the partial bond energy) E^(n) etc. 
asi7 

Eij = Wij/Z + E;(n)/n + E;(m)/m (1) 
where i, j = A or B atom of the general A^Bx-i alloy and w^ = if i = j. 
Ep(n) is the partial bond energy of the i-atom with n nearest neighbors etc. w^ 
is the interchange energy for two dissimilar atoms which may be estimated from 
the excess heat of mixing; and Z is the bulk coordination number. In order to find 
Eg(n) Strohl and King used an empirical formula 

E^(n) = a + b'n + cn (2) 
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where a , b and c are constants for the i-metal and nnay be obtained from the 
experimental values of the heat of sublimation, the energy of bulk vacancy formation 
and the surface energy for (100) surface. 

We use the a, b and c values evaluated by Rousset et al.^^ from a modified 
tight binding method (MTB). Here, the parameters are determined from the ex- 
perimental values of dimer energy, surface energy and the energy for monovacancy 
creation. In this MTB scheme the site energy for six-coordinated corner atom and 
the seven-coordinated edge atoms can be obtained from interpolation rather than 
extrapolation as was obtained in the work of Strohl and King. Once the site en- 
ergies are known, the pair bond energies can be easily evaluated. The second 
modification is made to consider the chemisorption. For the present study we have 
introduced necessary changes in the MC simulation to take into account the effect 
of chemisorption. Here we take an additional term 6{Ea-^b) in the configuration 
energy where 6 is the adsorbate coverage and Ea{b] is the chemisorption energy of 
the adsorbate on A{B) metal. Once the constants a , 6* and c are derived 
and the parameters w^/Z, 9 and (E^-Eg) are provided, the MC simulation is car- 
ried out in the usual manner until the equilibrium configuration energy is found. In 
essence, the atomic distribution results obtained in our nanostructure work would 
complement the segregation results obtained earlier by Rousset et al.^^ for the 
semi-infinite Pd-X alloy systems. In addition, the present results consider the role 
of various degrees of coverage of adsorbates on the atomic distribution. 

III. APPLICATIONS 

In view of the fact that all the constituent metals of the bimetallic systems con- 
sidered in this work have fee structure, the particles with number of atoms in the 
range of 200 - 2500 may be assumed to have the fee eubo-octahedral geometry. 
With this assumption we perform the MC simulation. The important input parame- 
ters used in the work are presented in Tables 1 and 2. In this simulation program the 
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important variable parameters are (a) the bimetallic system, (b) the bulk composi- 
tion, (c) the adsorbate coverage, and (d) the temperature. Since the composition 
variation with temperature is well-known, we have performed the simulation for 
a typical alloy equilibration temperature (700K). Though we have performed sim- 
ulation for several bulk compositions (20% Pd, 50% Pd, 66% Pd, 80% Pd), for 
convenience of comparison of several systems in presence of different adsorbates 
with varying coverage we would present results mostly for PdsoXso composition. 
The simulation results being very exhaustive cannot all be presented in this paper. 
We would be presenting some typical results of the systems. Readers interested in 
the details for other compositions, particle sizes, coverage or adsorbates may kindly 
contact the authors. 

A. Average surface Pd concentration 

In figure 1(a) we have plotted for the PdsoAgso nanosystem the fraction of 
surface sites occupied by Pd atoms for different particle size and H coverage. In 
Figs. 1(b), (c) and (d) similar results are plotted for the oxygen, CO and NO 
adsorption. Several characteristics are to be noted: 

(i) For clean Pd5oAg5o nanoparticles of all sizes Ag segregates to the surface sites. 

(ii) As the coverage increases the fraction of Pd surface atom increases leading to 

a reduction in Ag segregation and ultimately Pd atoms tend to segregate to 
the surface sites. 

This is true for all particle sizes and all the four adsorbates under consideration. 
In case of CO adsorbate Pd atoms seem to segregate to the surface sites even at 
coverage 0.25, while for other gases the segregation reversal takes place at coverage 
close to one. This behavior can be interpreted as a result of the higher heat of 
adsorption of adsorbates on Pd compared to that on Ag. Higher binding energy of 
adsorbates causes Pd atoms to move to the surface countering the normal segregat- 
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ing force. The difference in heat of CO adsorption on Pd and Ag being very large 
(f«28 kcals/mole), Pd atoms start segregating to the surface even at low coverage. 
It may be noticed that for clean and small coverage of adsorbates the Ag surface 
fraction increases as the particle size increases (smaller dispersion). This follows 
from the fact that the fraction of surface sites decreases with increase in particle 
size; and therefore, the Ag atoms, if available, can cover higher fraction of surface 
sites for clean and small coverage (please see Table 4 of Strohl and King^^). How- 
ever, as the coverage increases, more Pd atoms move to the surface due to stronger 
binding to the adsorbate; and thus, the Pd surface fraction increases with increase 
in particle size. 

Figs. 2(a)-(d) show the Pd surface fraction for the PdsoCuso nanoparticles in 
presence of the adsorbed gases. For this bimetallic system, Cu segregates to the 
surface for hydrogen coverage up to 0.6 for larger particles. For smaller particles 
Cu segregates for hydrogen coverage even up to 0.75. For still larger coverage 
Pd segregates to the surface. For oxygen as adsorbates on Pd5oCu5o system, Cu 
segregates for all coverages and all particle sizes studied in this work. Higher the 
coverage, higher is the extent of Cu segregation. This again follows from the heat 
of adsorption difference criterion mentioned before. Here binding energy of oxygen 
on copper is higher than that on Pd. For CO and NO adsorption on the PdsoCuso 
nanoparticles Cu segregates for clean and very low coverage. For CO and NO 
coverage >0.25 the Pd atoms start segregating to the surface. 

The difference in the heat of adsorption of hydrogen (oxygen) on pure Pd and Ni 
metals is very small. This causes only small changes in the surface Pd concentration 
with coverage for the PdsoNiso system (Figs. 3(a) and (b)). For all particle sizes 
and coverage up to one Ni segregates to the surface. For CO and NO adsorption 
on PdsoNiso, on the other hand, Pd atoms start segregating to the surface from 
coverage 0.2 (please see Figs. 3(c) and (d)). The segregation behavior of the 
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PdsoPtso system is very simple. The clean system shows Pd segregation. And for 
all particle sizes and all adsorbates Pd surface concentration slowly increases with 
increase in coverage (please see Figs. 4(a)-(d)). 

In order to compare the segregation behavior of the bimetallic PdsoXso nanosys- 
tems in presence of different adsorbates we plot in figure 5(a)-(d) the Pd surface 
fraction for a particular coverage, 9 — 0.5. It may be noticed that in presence of 
hydrogen or oxygen (Figs. 5(a) and (b) respectively) Pd segregates to the surface 
sites only for the Pd-Pt system. For the rest of the systems the other compo- 
nent segregates. In presence of CO, however, Pd segregates for all systems and 
for all particle sizes (Fig. 5(c)). In presence of NO, only for the Pd-Ag system 
Ag segregates, while for the rest of the systems Pd segregates to the surface sites 
(Fig.5(d)). The role of adsorbate coverage on the surface Pd concentration for the 
different bimetallic systems has been shown in Figs. 6(a)-(d) for the four different 
adsorbates. The results are for particles having 586 atoms (dispersion D = 0.457). 
The figures are self-explanatory. Except for hydrogen on Pd5oNi5o (Fig. 6(a)) and 
oxygen on PdsoNiso and PdsoCuso (Fig. 6 (b)) for all other systems Pd surface 
fraction increases with coverage. 

B. Site occupancy and surface coordination 

Besides the overall surface Pd concentration for Pd-X nanoparticles it is also very 
important to know the occupancy of various types of surface sites on the fee cubo- 
octahedron particle. On the particle surfaces there are 6-coordinated corner sites, 
7-coordinated edge sites, 8-coordinated fcc(100)-like face sites and 9-coordinated 
fcc(lll)-like face sites. In figures 7(a)-(d) we have shown how for CO adsorption 
on 586-atom nanoparticles of the four bimetallic systems the fraction of Pd atoms 
on the 6, 7, 8 and 9-coordinated surface sites changes with coverage. For CO on 
PdsoAgso system (Fig. 7(a)) at low coverage the 6 and 7 -coordinated sites are 
generally occupied by Ag atoms. As the coverage increases, gradually all types of 
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sites get increasingly occupied by Pd atoms. And at full monolayer coverage except 
the 8 -coordinated sites all other sites are occupied by Pd atoms. This is also the 
case for PdsoCuso (Fig. 7(b)) system. In case of PdsoNiso system (Fig. 7(c)) 
no particular type of surface sites is occupied fully by Pd or Ni atoms at any of 
the coverage from zero to one. In case of PdsoPtso system the occupancy of the 
9-coordinated sites changes very slowly from 0.72 to 0.78 as the coverage increases 
from zero to one. For NO adsorption the coverage-dependent site occupancy results 
for the four bimetallic systems are shown in figures 8(a) to (d). The results are self- 
explanatory. 

For studying the relative occupancy of various sites in the four bimetallic sys- 
tems we plot in figures 9(a)-(b) the site occupancy in the systems for a particular 
coverage (6) and a particular particle size (586-atom nanoparticles). The adsor- 
bates considered are CO and NO. For both the adsorbates, it may be noticed that 
at this 0.5 coverage the 6-coordinated sites are almost fully occupied by Pd in case 
of Pd5oPt5o system. On the other hand, for the Pd5oAg5o and PdsoCuso systems 
the 6 -coordinated sites are generally occupied by Ag and Cu atoms respectively. 

C. Size dependence of surface site occupancy 

Particle size, quite often, plays an important role in catalysis. In bimetallic 
systems, this has an additional significance. Because, different metallic atoms may 
preferentially occupy different type of sites. In order to study this aspect we show 
in figure 10(a)-(d) how for a CO coverage of 0.5 the surface sites of a 586-atom 
nanoparticle are preferentially occupied by one type of atom or the other. What we 
find is that there is no general rule for occupancy of various sites as the particle 
size increases as there are crossovers of curves for different particle sizes in each of 
the four systems. In case of Pd5oAg5o system (Fig. 10(a)) the occupancy of 8 and 
9-coordinated sites by Pd atoms does not change much with particle size. But there 
are fluctuations in the occupancy of 6 and 7-coordinated sites by Pd atoms as the 
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particle size increases. Similarly, in case of PdsoCuso and PdsoNiso system there are 
fluctuations in occupancy of all types of surface sites as the particle size increases 
(Figs. 10(b) and (c)). In case of PdsoPtso we find some systematic variation in the 
occupancy of various surface sites. As the number of atoms in the particle increases 
from 201 to 2406, the occupancy of 9-coordinated sites by Pd atoms increases 
significantly. This variation is distinctively different from other bimetallic systems. 

The net result of the model presented in this work is that we have simulated 
the probable surface composition and the atomic distribution for the nanoparticles 
of some Pd-based bimetallics in presence of various coverages of simple gases like 
hydrogen, oxygen, CO and NO. However, the experimental results are scarce, par- 
ticularly for adsorbed systems. In the zero-coverage limit, i.e. for clean bimetallic 
nanoparticles experimental results exist for the PdPt and PdCu systems. The model 
shows Pd segregation for the PdPt system and Cu segregation for the Pd-Cu system 
in agreement with the experimental results. ^^'^^ But the coverage-dependent sur- 
face composition and site occupancy are yet to be experimentally realized. However, 
in the absence of experimental results, the present model can serve some purpose 
in estimating the probable surface composition for the nanoparticles under various 
coverage of gases. Besides, the knowledge of the site occupancy by A type atom or 
by B type atom of A-B bimetallic atoms may sometime help to identify the active 
sites on a catalyst for a particular catalytic reaction. For example, it has been found 
that the rate of 1,3-butadiene hydrogenation reaction on a Pd5Ni95 system is larger 
than that on pure metals. This may be due to the ensemble effect whereby a 
particular form of ensemble of Pd and Ni atoms form the active site. Such behavior 
may be explained from a knowledge of the site occupancy as studied in this work.^'' 

IV. CONCLUSIONS 

Monte-carlo simulations are carried out to find out the nature of occupancy 
of surface sites of Pd-based bimetallic nanoparticles. We upgraded the earlier MC 
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simulation scheme of Strohl and King^^ with input data obtained from improved 
cohesive energy calculation.^^ 

(a) It has shown how the overall surface composition varies with particle 
size, the nature of adsorbed gas and its coverage . 

(b) It has also been shown how for a particular adsorbed gas the surface 
composition varies from one bimetallic system to the other. 

(c) We have also shown how different types of surface sites of the bimetallic 
systems are occupied and how the occupancy of these sites varies with 
gas coverage. 

(d) We have discussed briefly the importance of such studies in understand- 
ing the role of active sites in catalysis. 

The validity of the model would be justified after more and more experimental 
data and particularly for gas-covered surfaces, are available. 
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FIGURE CAPTIONS 



Fig. 1. MC simulated surface fraction of Pd atoms on bimetallic PdsoAgso 
nanoparticles as a function of dispersion and coverage of various adsor- 
bates: (a) Adsorbate: hydrogen; (b) Adsorbate: oxygen; (c) Adsorbate: 
CO; (d) Adsorbate: NO. 

Fig. 2. IVIC simulated surface fraction of Pd atoms on bimetallic PdsoCuso 
nanoparticles as a function of dispersion and coverage of various adsor- 
bates: (a) Adsorbate: hydrogen; (b) Adsorbate: oxygen; (c) Adsorbate: 
CO; (d) Adsorbate: NO. 

Fig. 3. MC simulated surface fraction of Pd atoms on bimetallic PdsoNiso 
nanoparticles as a function of dispersion and coverage of various adsor- 
bates: (a) Adsorbate: hydrogen; (b) Adsorbate: oxygen; (c) Adsorbate: 
CO; (d) Adsorbate: NO. 

Fig. 4. MC simulated surface fraction of Pd atoms on bimetallic PdsoPtso 
nanoparticles as a function of dispersion and coverage of various adsor- 
bates: (a) Adsorbate: hydrogen; (b) Adsorbate: oxygen; (c) Adsorbate: 
CO; (d) Adsorbate: NO. 

Fig. 5. MC simulated surface fraction of Pd atoms on bimetallic PdsoXso 
nanoparticles as a function of dispersion and 0.5 coverage of various 
adsorbates. (X = Ag, Cu, Ni, Pt). (a) Adsorbate: hydrogen; (b) 
Adsorbate: oxygen; (c) Adsorbate: CO; (d) Adsorbate: NO. 

Fig. 6. MC simulated surface fraction of Pd atoms on 586-atom bimetallic 
PdsoXso nanoparticles as a function of coverage of various adsorbates. 
(X = Ag, Cu, Ni, Pt). (a) Adsorbate: hydrogen; (b) Adsorbate: oxygen; 
(c) Adsorbate: CO; (d) Adsorbate: NO. 
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Fig. 7. MC simulated surface fraction of Pd atoms at N-coordinated surface 
sites (N= 6-9) of a 586-atom PdsoXso nanoparticle in presence of var- 
ious coverage of CO. (a) X= Ag; (b) X= Cu; (c) X=Ni; (d) X= Pt. 

Fig. 8. MC simulated surface fraction of Pd atoms at N-coordinated surface 
sites (N= 6-9) of a 586-atom PdsoXso nanoparticle in presence of var- 
ious coverage of NO. (a) X= Ag; (b) X= Cu; (c) X=Ni; (d) X= Pt. 

Fig. 9. MC simulated surface fraction of Pd atoms at N-coordinated surface 
sites (N= 6-9) of a 586-atom PdsoXso nanoparticle (X — Ag, Cu, Ni, 
Pt) for 0.5 coverage of adsrbates. (a) Adsorbate: CO; (b) Adsor- 
bate:NO. 

Fig. 10. MC simulated surface fraction of Pd atoms at N-coordinated surface 
sites (N= 6-9) of PdsoXso nanoparticle of different sizes in presence of 
0.5 coverage of CO. (a) X= Ag; (b) X= Cu; (c) X= Ni; (d) X= Pt. 
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Table 1 

Input energy parameters for Pd-X bimetallic particles 



system 




Eh 


Eo 


Eco 


Ejvo 


Pd-X 


(ev) 


(kcals/mol) 


(kcals/mol) 


(kcals/mol) 


(kcals/mol) 


Cu 


-0.0197 


-56" 


-103° 


-12° 


-14' 


Ag 


-0.00957 


-56° 


-80° 


-6° 


-25° 


Ni 


-0.0095 


-63° 


-90' 


-27° 


-25° 


Pt 


-0.00396 


-61" 


-85'' 


-32° 


-27° 


Pd 





-62° 


-87° 


-34° 


-31° 



Ref. 22, Ref. 23, Ref. 24, Ref. 25, Ref. 15 



Table 2 

a, b, and c parameters for different metals 





P(l 


Cu 




Ni 


Pt 


a 


-0.17702 


-0.34040 


-0.25866 


-0.42578 


-0.42874 


b 


- 0.04842 


-0.01177 


-0.013283 


-0.01447 


-0.04750 


c 


0.00299 


0.00129 


0.00119 


0.00160 


0.00355 
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